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We study and characterize a direct route to high- dimensional chaos (i.e. not implying an interme- 
diate low-dimensional attractor) of a system composed out of three coupled Lorenz oscillators. A 
geometric analysis of this medium-dimensional dynamical system is carried out through a variety of 
numerical quantitative and qualitative techniques, that ultimately lead to the reconstruction of the 
route. The main finding is that the transition is organized by a heteroclinic explosion. The observed 
scenario resembles the classical route to chaos via homoclinic explosion of the Lorenz model. 
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Low-dimensional dissipative continuous sys- 
tems, namely described by a d=3 phase space, 
may exhibit chaotic behavior, that may arise 
through a number of well characterized scenarios, 
including the route from quasiperiodicity to chaos 
and routes involving global bifurcations, like the 
route to chaos exhibited by the Lorenz system. 
Things become more complicated when the di- 
mensionality of the phase space is increased. In- 
terestingly, the limit in which one has an infinite- 
dimensional phase space, and the system exhibits 
the so-called spatio-temporal chaotic behavior, is 
well characterized through the ergodic theory of 
dynamical systems, through quantifiers like di- 
mension, entropy, and Lyapunov exponents (a 
Lyapunov density in continuous systems). Less 
explored is the case in which one has a phase 
space that has an intermediate dimension be- 
tween these two limits. In this case, one may 
have high-dimensional chaotic behavior, and an 
obvious possible route is a cascade in which the 
dimensionality of the high-dimensional attrac- 
tor increases sequentially. Still, other possibili- 
ties may take place, like the scenario discussed 
in the present work, in which quasiperiodic at- 
tractors participate, though in a manner differ- 
ent to the classical quasiperiodic route to chaos. 
Instead, a high-dimensional chaotic attractor is 
created in a global bifurcation. Interestingly, 
the route exhibits a regime in which a three- 
dimensional quasiperiodic attractor is stable, in 
apparent contradiction to common wisdom inter- 
pretation of Ruelle-Takens-Newhouse Theorem. 
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I. INTRODUCTION 

Nonlinear dynamics was considerably boosted by the 
recognition of chaos as an ubiquitous phenomenon out 
of the linear realm. Chaos theory succeeded to iden- 
tify the universal routes through which regular motion 
may become chaotic. In dissipative dynamical systems, 
chaos appears through a few well characterized routes (or 
scenarios) [H like the following: (a) the period-doubHng 
cascade [2'|; (b) the intermittcncy route [Sll; (c) the route 
involving quasiperiodic tori [4,, ^] ; (d) the crisis route [1) ; 
(see also 0,11] for a survey). Another possibility is that 
chaos sets in through a global connection to a fixed point, 
as is the case [9] of the Lorenz system, and of Shilnikov 
chaos US [Hi. 

More recently, some studies have been devoted to the 
study of high-dimensional chaos, that here we shall de- 
fine as chaotic motion inside an attractor of typical di- 
mension d > 3 (as a typical one we may take the infor- 
mation dimension Di). An obvious possibility to transit 
to high-dimensional chaos is to take as starting point a 
low-dimensional chaotic attractor. This possibility has 
been often considered in the context of desynchroniza- 
tion between coupled chaotic oscillators. A possible sce- 
nario, and probably the most studied one [H, [l^ in- 
volves generating a hyperchaotic attractor — an attractor 
with two or more positive Lyapunov Exponents (LEs) — 
from a low-dimensional chaotic attractor. Another class 
of high-dimensional chaotic attractors exhibits two (at 
least very approximately) null Lyapunov exponents. As 
shown in T4, 15, 16] one may encounter this situation 
in rings of asymmetrically (e.g. unidirectionally) coupled 
chaotic systems in which a symmetric Hopf bifurcation 
renders unstable the synchronized state [ij, [l^ . 

Less obvious is the possibility of a direct transi- 
tion to high-dimensional chaos without an intermedi- 
ate low-dimensional chaotic attractor. Restricting to 
autonomous ordinary differential equations, the transi- 
tion to high-dimensional chaos has been found to be 
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associated to quasiperiodicity. Thus, the works by 
Feudel et al. il9j and Yang [l^ report a transition 
from two- and/or three- frequency quasiperiodicity to 
high-dimensional chaos. A more geometrical view is 
in the work by Moon that describes a route (to 

high-dimensional chaos) from a two-dimensional torus, 
through a global bifurcation that comprises a double 
homoclinic connection to a limit cycle (this structure 
amounts to adding one dimension to each building block 
of the homoclinic route to chaos in the Lorenz model). 
Recently [2l|, we reported a route to chaos whose most 
novel aspect is that it implies the sudden creation (i.e. 
with no mediating low-dimensional chaos) of a high- 
dimensional chaotic attractor with Di > 4. In summary, 
in this route first a high-dimensional chaotic (nonattract- 
ing) set is created in a global bifurcation. Then, this set 
is rendered stable, becoming a high-dimensional chaotic 
attractor (56j . in a boundary crisis. 

The purpose of this paper is to explain in more detail 
our findings in Ref. [2l|. In Sec.|TT]the system is described 
and an overall picture of the route to high-dimensional 
chaos is presented. Then, Sec. lIIII discusses at some detail 
the appearance and the robustness of the two- and three- 
frequency quasiperiodic attractors found in the system. 
Later, Sec. lIVI presents the numerical evidences that have 
been used to understand the complex route to chaotic 
behavior presented by the system. Section |V] presents a 
characterization of the route to high-dimensional chaos 
through a return map similar to that of the Lorenz sys- 
tem. The goal of Section IVII is to present and discuss at 
some length the route to chaos exhibited by the system. 
Finally, Sees. IVIJ and rvTIIl present the final remarks on 
this work and the conclusions, respectively. 



II. SYSTEM AND OVERALL PICTURE 

The system studied in the present work is a 9- 
dimensional dynamical system considered here is formed 
by three Lorenz oscillators coupled according to the 
partial replacement coupling method of Ref. [IH . The os- 
cillators are unidirectionally coupled in a ring geometry 
such that. 



spatial modes [17|,l27|. These modes are defined as, 
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~ Xj-i for J ^ 1, define the coupling, and 
xi = the ring geometry of the system. In our study 
two of the parameters are fixed as in [U [ISl)- c = 20, 
6 = 3, while 29 < i? < 40. The study of this system has 
been suggested by the results of the experimental study 
of rotating waves for three coupled Lorenz oscillators cor- 
responding to these parameters [H, . 

A useful representation in the study and characteriza- 
tion of discrete rotating waves are the (discrete) Fourier 
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where = 3 (as already indicated) and i is the imaginary 
unit. In terms of these modes (Xg G R"^,Xi G C'^), the 
evolution equations are. 



Ao = a(yo-Ao) 

Ya = R Ao —Yq— Aq — Xi Zl — XI Z\ 

Zo = Ao Yo - ^0 + Al Yi* + XI Yi 

Al = fT(Yi - Al) 

Yi = — Yi — ^0 -^1 — ^1 -^0 — ^1 ■^1 

Zl = Ao Yi + Al Yo - 6 Zl + Ai* Y( 
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with R — Rcxp{2Tri/3), and where Ai denotes the com- 
plex conjugate of Ai. 

The qualitative description of the behaviors exhibited 
by the system in this range of parameters is as follows. 
The three Lorenz oscillators exhibit chaotic synchroniza- 
tion for R < Rsc ~ 32.82. At Rgc the system exhibits a 
Hopf bifurcation directly from a chaotic state, yielding a 
behavior, that was first found in rings of unidirectionally 
coupled Chua's oscillators and called a Chaotic Rotating 
Wave (CRW) [H]. The Hopf bifurcation exhibited by the 
system is called symmetric (28l . [IP] , as it originates from 
the cyclic (Z3) symmetry of the ring. Close to the bifur- 
cation the CRW is the combination of the chaotic dynam- 
ics in a Lorenz attractor and a superimposed oscillation 
created by the Hopf bifurcation, that occurs in the sub- 
space transverse to the synchronization manifold (fc = 1 
mode). The oscillation created in the symmetric Hopf 
bifurcation is characterized by a phase difference of 27r/3 
{2tt/N in the general A-oscillator case) between neighbor 
oscillators, and this behavior is the discrete analog of a 
traveling rotating wave [l^l . This picture is valid close to 
onset, R > Rgc, and when increasing R one observes that 
this behavior changes. For R > 35.26 the behavior of the 
system becomes periodic, with a waveform characteristic 
of the fc = 1 Fourier mode. Actually, two mirror limit- 
cycle attractors, called Periodic Rotating Wave (PRW) 
in Ref. [13, [Mi, are found in the range R G [35.26, 39.25] 
(approximately). At i? = RpUch ~ 39.25 both solutions 
merge giving rise, through a pitchfork bifurcation, to a 
centered stable symmetric periodic behavior. 

We shall not describe here the transition from syn- 
chronous chaos to CRW; instead, we shall focus on the 
transitions from PRW to CRW, i.e. we go 'from order 
to chaos' decreasing R. The reader should notice that 
by CRW we refer to a high-dimensional chaotic attractor 
characterized by a oscillation with a 27r/3 phase shift be- 
tween neighboring oscillators superimposed to an under- 
lying chaotic behavior. The temporal series for the dif- 
ferent behaviors studied in this work are shown in Fig. [1] 
Note that the third frequency, that appears as the T'^ is 
born, manifests as a very slow modulation on the size 
of the former (the time scale has been broadened in 
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FIG. 1: Time series of the system Eq. (HJ: (a) R = 35.5, Pe- 
riodic Rotating Wave; (b) R = 35.2, T^; (b) R = 35.095, T^; 
(d) R = 35.093, Chaotic Rotating Wave. Note the different 
time scale for each panel. 



Fig. [Ijc) in order to allow the observation of this slow 
scale). 

A view of the parameter region that will be considered 
in this work can be found in Fig. [21 in which the Lya- 
punov spectrum corresponding to the four largest Lya- 
punov exponents is presented. From right to left, the two 
symmetry-related limit-cycle attractors exhibit a (super- 
critical) Hopf bifurcation when R is decreased, namely 
at i? — w 35.26, yielding two symmetry-related two- 
frequency quasiperiodic attractors (with two null LEs). 
Lowering R further the system exhibits another (super- 
critical) Hopf bifurcation, at i? = Rh2 ~ 35.0955, that 
yields two mirror three-frequency quasiperiodic attrac- 
tors (three vanishing LEs). Lowering R furthermore 
the system exhibits a boundary crisis, at R — Rtc ~ 



35.09384, in which the chaotic attractor is born (or de- 
stroyed seen from the opposite side). The chaotic at- 
tractor is characterized by two vanishing (at least, very 
approximately) Lyapunov exponents, and a single pos- 
itive LE that is larger than the absolute value of the 
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FIG. 2: The four largest Lyapunov exponents (Ai=i,...,4) as 
a function of the parameter R. Four regions can be distin- 
guished according to their different Lyapunov spectra. Note 
that for the high- dimensional chaos we obtain A2 = A3 = 
and Ai > IA4I which implies, according to the Kaplan- Yorke 
conjecture, an information dimension Di > 4. The dash- 
dotted line indicates the value of R where the chaotic attractor 
is born. 



fourth LE. This implies, according to the Kaplan- Yorke 
conjecture 57|, that the information dimension Di > 4. 
As we shall see below (see Sec. IIVPP this dimension is 
genuine. 



The presence of 2-D and 3-D quasiperiodic attrac- 
tors may lead to think that chaos appears through a 
quasiperiodicity transition to chaos (see Sec. IIII|) . How- 
ever, it will be shown below that the chaotic attractor is 
created at a boundary crisis, and coexists with the two 
attractors until the latter are destroyed as each of them 
collides with a twin unstable T'^, at i? = Rgn ~ 35.09367. 
We shall show (from Sec. lIV[) that, indeed, the system ex- 
hibits a global bifurcation — analogous to what happens 
for the Lorenz system Q — that implies, in first approx- 
imation, the sudden creation of an infinite number of 
unstable 3-D tori. 

A schematic diagram of the whole set of bifurcations 
linking the PRW and synchronous chaos is shown in 
Fig. [3] As mentioned above, in the interval of R where 
the CRW is found, the shape of the attractor changes. 
Anyway, in this paper we are not interested in the tran- 
sitions between different types of chaotic rotating wave 
(see [15 for such a study). 



III. TRANSITION TO QUASIPERIODIC replaced after a sequence of bifurcations by two mirror 

BEHAVIOR 

As described above, the (symmetric) periodic rotat- 
ing wave, that is the starting point of our analysis, is 
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FIG. 3: Diagram representing schematically the transitions from synchronous chaos (left) to a PRW (right). 



three- frequency quasiperiodic attractors, T''. The in- 
volved transitions are a pitchfork bifurcation, in which 
two (asymmetric) periodic rotating waves are born, and 
then two consecutive Hopf bifurcations. A small note 
is in place about the stability of this type of unusual 
attractors as it has been stated that T'^ attractors are in- 
trinsically fragile according to Newhouse-Ruelle-Takens 
(NRT) Theorem [i]. 

First of all, one can present numerical evidence to show 
that, indeed, three- frequency quasiperiodic attractors are 
observed robustly in a finite parameter range. In partic- 
ular the inset of Fig.[2]shows results corresponding to the 
calculation of the Lyapunov spectrum in the T'^ param- 
eter range. As there are long chaotic transients in this 
region, the calculations were performed using long tran- 
sients (of the order of 10^ t.u.), while another 10^ t.u. 
were used to calculate the spectrum. 

On the other hand, and as already reported in Ref . [13] , 
it can be numerically checked that the system exhibits 
three-frequency quasiperiodic behavior by looking at the 
Poincare cross section of a T'^ torus in phase space. 
Namely, in Fig. 2 of Ref. [24| the Poincare sections of a 
and of a T"^ torus, respectively, are presented. While the 
former is a circle, the latter is clearly a two-dimensional, 
T^, torus. In the literature Ref. 16] reports the existence 
of three-frequency quasiperiodicity in a ring of coupled 
Lorenz oscillators (in a different parameter regime) , while 
a case of a third very low frequency associated to a T"^ 
torus has been observed in Refs. [33. Issl. [sl] . 

In trying to justify the robust existence of T'^ dynam- 
ics in the light of the NRT Theorem, a first observa- 
tion is that in our simulations we have not detected reso- 
nances in the regime. This may occur due to the fact 
that along i?, the winding number does not crosses hard 
resonances, the smallest denominator is g = 11 (corre- 
sponding to a winding number of 4/11), so the resonance 
horns may be quite small. But also it is important to no- 
tice that the absence of resonances in rotating waves is a 
characteristic feature of systems with rotational, 50(2), 
symmetry 35] . For instance in a homogeneous excitable 
medium, the compound rotation of a spiral wave lacks of 
frequency lockings due to this symmetry [36|. Therefore 
the stability of the T'^ attractors would be subsequent to 
the robustness of the quasiperiodic dynamics on the two- 
tori. In a sense our three-frequency quasiperiodic attrac- 



tors would be as robust as two-frequency quasiperiodic 
attractors. It is clear that in our system this rotational 
symmetry is only approximate, but one may consider the 
discreteness effect as a small perturbation. Also the third 
frequency is very small what precludes strong resonances 
with small denominators: thus, for R = 35.095 [Fig.[IJc)] 
we have h = 5.3 . . ., /s = 1.87 . . ., and h = 0.016 . . .. 
These arguments support the existence of a range of 
three-frequency quasiperiodic behavior. 



IV. NUMERICAL EVIDENCES OF THE ROUTE 
TO CHAOS EXHIBITED BY THE SYSTEM 

In this section we are going to analyze the route 
through which a chaotic attractor is born in this sys- 
tem. First of all, we shall reduce the dimensionality of 
the problem by eliminating the fast frequency, namely the 
one involved in the phase shift by 27r/3 in neighbor oscil- 
lators, as it leads to a conserved quantity (this time lag) 
and consequently a vanishing (or almost vanishing) Lya- 
punov exponent. In the mode representation of Eq. ([3]) 
this amounts to perform a cut through the Poincare sec- 
tion Im{Xi) = 0, Jm(Xi) > 0. So, in visualizing objects 
cycles will become fixed points, T^-tori cycles and T^-tori 
will become T^-tori, although we shall refer to these ob- 
jects corresponding to the complete phase space, and not 
to the result of the dimensionality reduction achieved via 
sectioning. 



A. Coexistence between 3D-torus and CRW 

The first important remark about our system of three 
coupled Lorenz oscillators is that the two T'^ attrac- 
tors are not directly involved in the birth of the high- 
dimensional chaotic attractor. Continuation of the T'^ 
attractors shows that they exist above Rsn = 35.09367. 
This implies that the system exhibits multistability be- 
tween the high-dimensional chaotic and the T^^ attractors 
in the range Rsn < R < Rbc, with the latter attractors 
having a small basin of attraction. This remark is impor- 
tant, because it implies that the high-dimensional chaotic 
attractor is not created through some route involving the 
T"^ attractors. 
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FIG. 4: Blowout of the four largest Lyapunov exponents for 
the attractors in the region in which they coexist with the 
chaotic attractor. The fourth Lyapunov exponent approaches 
zero in a square-root fashion, as expected for a saddle-node 
bifurcation. 

A detailed view of the behavior of the four largest Lya- 
punov exponents for the T'^ attractors is represented in 
Fig. m It can be seen that close to Rsn the fourth Lya- 
punov exponent exhibits a square-root profile as expected 
for a saddle-node bifurcation. This suggests that each 
attractor is approaching an unstable T^. Also one can 
appreciate quite clear lockings [5^] where the third and 
fourth Lyapunov exponents become equal. 

As R approaches Rsn the existence of a smooth invari- 
ant 3D-torus cannot be guaranteed, because the strength 
of the normal contraction can be of the same order than 
the rate of attraction of the 2D-torus arising in each lock- 
ing [33|. Thus, in Fig. [5] a possible scenario for the in- 
teraction of the tori, that would explain the lockings ob- 
served in Fig. m is presented. Because of the extremely 
small transversal stability of the T'^ (represented by the 
fourth exponent), for R > Rsm we conjecture that the 
transversal direction 'slaves' the tangential one (repre- 
sented by the third exponent). Then, we observe two 
identical non-vanishing exponents that may indicate the 
existence of a stable focus-type on the surface of the 

. The T'^ continues to exist but is non-differentiable at 
the stable located on its (hyper)surface. 

The participation of the unstable 3D-torus in the final 
annihilation of the stable one is not trivial (the existence 
of the unstable is supported by another numerical evi- 
dence, discussed in Sec. lIVE)) . but we want to recall that 
saddle-node bifurcations of T-^ tori have been reported 
[38] (called bubbles in this reference), and due to the 
special character of the T'^ tori that we have here (with 
one oscillation that, in some sense, is orthogonal to the 
others and also almost conserved, cf. Sec. IIIip . and given 
the numerical evidence reported here, we think that we 
have found a saddle-node bifurcation of T'^ tori. 



FIG. 5: Schematic of a 'slave locking'. A portion of a section 
of the tori (stable and unstable) are shown with bold lines. 
Tori are not differentiable at the foci. 



B. Route to chaos in the Lorenz model through 
homoclinic explosion and boundary crisis 

As the route studied in this work presents some sim- 
ilarities with the classical route to chaos for an isolated 
Lorenz system, we are going to draw some useful analo- 
gies with that system, although it is warned that the 
analogy is not complete (otherwise, we would not have a 
new route to chaos). As it can be found in textbooks [1], 
the Lorenz system exhibits several routes to chaos. The 
classical one is the route to chaos through a double homo- 
clinic connection of the saddle point at the origin. Fixing 
a and b, the double homoclinic orbit occurs at a particu- 
lar value R = Rho m [E^ • It has the shape of the butterfly 
(taking into account how the unstable directions reenter 
the saddle point at the origin, being this dictated by the 
reflection symmetry of the system) , and this gives birth 
(lol | to a chaotic set for R > Rhom (a homoclinic ex- 
plosion 0, H^). The closure of this set is formed by the 
infinite number of unstable periodic orbits (UPOs) that 
can be classified according to their symbolic sequences of 
turns around the right (R) and the left (L) fixed points 
(C+ and C_) [9j|. The appearance of these UPOs reflects 
the dramatic change undergone by the stable manifold of 
the flxed point that allows initial conditions at one side of 
the phase space jump to the other side (before falling to 
one of the two symmetry related stable fixed points C+ 
and C_), being these jumps impossible for R < Rhom- 
For a larger value of R, R = Rbc, the chaotic set be- 
comes stable in a boundary crisis, that occurs precisely 
when the two shortest symmetry- related length- 1 unsta- 
ble periodic orbits (R and L) located at each lobe, and 
that at i? = Rhom coincide with the two homoclinic 
orbits, shrink, such that the chaotic set has a tangency 
with these two orbits At this value of R there ex- 
ists a double heteroclinic connection between the equi- 
librium at the origin and the mentioned length-1 UPOs. 
For Rbc < R < Rh these two orbits (and their respec- 
tive tubular stable manifolds) form the basin boundaries 
between the chaotic attractor and the two stable asym- 
metric fixed points (C±), and, consequently, the system 
exhibits multistability. These two fixed points loose their 
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FIG. 6: Numerical experiment showing the time series of the 
coordinate Xo at the Poincare section for trajectories starting 
in an initial condition in the symmetric unstable PRW: (a) for 
R = 35.12, a condition before the explosion (b) for R — 35.10 
just past the explosion. The Xo coordinate of one of the two 
asymmetric PRWs is denoted with a bold dotted line. 



stability in a a subcritical Hopf bifurcation, that occurs 
when the two mentioned length- 1 orbits shrink to a point 
coinciding with (C+, C_). 



C. Heteroclinic explosion 

In our system one can find a value of i? = Rexpi ~ 35.11 
that defines a clearcut transition in the way in which 
transients approach the attractors (that for R ^ Rexpi 
are two T^-tori). For R > Rexpi the basin of attraction 
of each is quite simple. But below Rexpi trajectories 
may tend asymptotically to one of the after visiting 
the neighborhood of the other torus. 

Following the analogy with the Lorenz system, we con- 
jecture that precisely at i? = Rexpi a global bifurcation 
occurs, and past this value (i? < Rexpi) an infinite num- 
ber of unstable objects are created. To check this we 
stabilized the symmetric PRW (an unstable fixed point 
in the Poincare map) by a Newton-Raphson method and 
observed the fate of the trajectories starting from (ap- 
proximately) that solution. In Fig.lHlthe evolution of the 
trajectories for two values of R above and below Rexpi is 
shown. After approaching one the asymmetric PRWs (its 
location is marked with a bold dotted line) the trajectory 
jumps or not to the other side. The result obtained for 
Rexpi is the same if one takes as starting point one of the 
asymmetric PRWs. 



D. Four-dimensional branched manifold 

The third important numerical remark refers to the 
Lyapunov spectrum (Table I) of the different attractors 



in the multistability region (of course, and due to sym- 
metry, the Lyapunov spectra of the two T'^ attractors are 
identical). As one can see, fifth to ninth LEs are quite 
similar, indicating that both attractor types are approx- 
imately embedded in the same four-dimensional space, a 
subset of the total phase space. 
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TABLE I: The five smallest Lyapunov exponents for the 
two attractors (CA: chaotic, and T'^: three- frequency 
quasiperiodic) coexisting at i? = 35.0938. 

Thus, we postulate that the dynamics can be simplified 
to the study of a four-dimensional branched manifold. A 
theoretical justification of this statement is impossible, 
but one can argue along the lines of (a higher-dimensional 
generalization of) the Birman- Williams Theorem [iO, |4l| 
(this theorem has proven to be very useful, even if it 
requires the strange attractors to be (uniformly) hyper- 
bolic what is usually not fulfilled; notice, however, some 
recent rigorous advances [42j on the theory of singular 
hyperbolic attractors, that encompass the Lorenz equa- 
tions f2l]. 

In our case we can argue that the high-dimensional 
chaotic attractor has, according to Kaplan- Yorke conjec- 
ture [43], a value of the information dimension Di > 4 
(cf. Fig. [2). We have also measured the correlation di- 
mension D2 (through Grassberger-Procaccia algorithm) 
obtaining also a value close to 4, namely D2 = 3.96 ± 
0.05), (recaU that Di>D2). 

A more correct picture of the attractor amounts to con- 
sider that this 4-dimensional manifold is actually com- 
posed of many thin leaves. Of course this reduced 4- 
dimcnsional picture can be only considered to be a more 
or less faithful representation of the system. In the 
epochs in which a trajectory jumps to the other side 
(subspace), which means reinjections through an extra 
dimension, it is when the existence of branching is needed 
in order to the trajectory not to intersect itself. In our 
case the 'tear point' is the symmetric PRW. This is in 
complete analogy with what happens with the Lorenz 
system, where the attractor can be understood as a tem- 
plate composed of a branched two-dimensional manifold 
with a tear point at the origin (cf. Fig. 20(b) in [i3|)- Ro- 
tations around one of the lobes are roughly planar, but 
reinjections between the two (planar) lobes, that form 
themselves an angle, involve the third dimension. 

An important remark concerning this 4-dimensional 
picture is that, in this space, a T"^ -torus divides a 4- 
dimensional manifold in two regions (just the same as 
a cycle divides a surface). Then in the regime with coex- 
istence, between chaos and three-frequency quasiperiod- 
icity, it is the pair of (conjectured) unstable 3D-tori what 
defines the basin boundary of the chaotic attractor. In 
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FIG. 7: Log-log representation of the average chaotic tran- 
sient as a function of the distance to the critical point Rtc- 
Each point is an average over 100 realizations. Data fit to a 
straight line of slope 7 — —1.53 ± 0.06. 
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FIG. 8: Numerical experiment showing a trajectory (only 
points at the Poincare section are shown) starting in an ini- 
tial condition at one of the asymmetric unstable PRWs for 
R — 35.09389 > Rbc- The unstable 3D-torus, larger, is seen 
initially, but finally (part (b)) the system decays to the smaller 
stable 3D-torus. Note that the minimum width along time is 
larger for the stable 3D-torus as expected for a torus located 
inside another (see the discussion about the branched mani- 
fold). 



the Lorenz system the length- 1 unstable periodic orbits 
divide the 2£'-manifold in three regions; and also they 
acts as the basin boundaries, when the chaotic and the 
fixed point attractors coexist. In the transient chaos re- 
gion, trajectories spiral in one lobe away from the fixed 
point (say C+) because of the repulsive effect of the UPO 
surrounding that fixed point. After some turns the tra- 
jectory jumps to the other lobe by using the third dimen- 
sion (i.e. by virtue of the branching). If the trajectory 
approaches close enough to C_ surpassing the 'barrier' 
constituted by the length- 1 UPO, it is captured by this 
fixed point and no more jumps occur (i.e. the chaotic 
transient finishes). Analogously, the unstable T'^-torus 
acts as a dividing hypersurface for trajectories in the 4- 
dimensional branched manifold, as discussed above. 



those R values, the chaotic transient is not observed, as 
may be seen in Fig. [5] (the PRW is a fixed point in the 
figure, as it has been stroboscopically cut through the 
Poincare hyperplane Im{Xi) = 0). For _R just above R^c 
(cf. Fig. [8|) the trajectory seems to approach a T'^ before 
'falling' to the T'^ attractor. Thus, we conjecture that is 
the unstable which constitutes the basin boundary of 
the chaotic attractor and the object involved in a global 
connection that marks the birth (or the death, depending 
on the viewpoint) of the chaotic attractor. Notice also, 
that the stable and the unstable T'^-tori have quite similar 
sizes. Ultimately, at Rsm the multistability region has 
an end when the two T'^-tori annihilate in a saddle-node 
bifurcation (actually there are two mirror bifurcations, 
as previously discussed in Sec. IIV A|) . 



E. Boundary crisis and power law of chaotic 
transients 

The fourth remark regards numerical studies for R > 
Rbc- We have measured the average time of the chaotic 
transients ((r)) when approaching R = R^c- We ob- 
serve that these transients diverge satisfying a power law 
((r) cx {R — Rbc)''), with 7 — — 1.53±0.06, for an asymp- 
totic value R = Rbc = 35.093838, as expected for a 
boundary crisis [ill, see Fig. Fffsot . 

Inspired in the boundary crisis occurring in the Lorenz 
system [H, we took one initial condition in one 
of the asymmetric (already unstable) Periodic Rotating 
Waves for different values of R. When R is above Rbc 
(but below Rexpi), the system goes away of the limit cycle 
(due to numerical roundoff errors), and jumps to the 
(or T'^) attractor located at the other side of phase space. 
Contrary to what happens to typical initial conditions for 



V. DESCRIPTION IN TERMS OF A RETURN 
MAP 

Lorenz [23| described a nice technique for reducing the 
complexity of the solutions of the Lorenz equations. By 
recording the successive peaks of the variable z{t), he 
reduced the dynamics of the Lorenz system to a one di- 
mensional map. Denoting the nth maximum of z(t) by 
Mn, he plotted successive pairs (M„, Mn+i), finding that 
points lay (very approximately) along a A-shaped curve. 
In this way, the dynamics is reduced to the "Lorenz 
map": M„+i = A(M„). 

In our case we may try to reduce the dynamics by 
means of a return map. The fast dynamics concerning 
the k ~ 1 mode is approximately filtered when consider- 
ing the fc = mode. Then, the T'^ attractor is seen in 
the k = framework as a plus some small residual 
oscillating component. A return map of the variable Zq 
would reduce the dimension of the attractors by one, and 
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FIG. 9: Return map of the maxima of the variable Zo satis- 
fying to be larger than their adjacent maxima (see text). The 
left panel shows a regime of transient chaos; after a transient 
the trajectory decays to the stable T'^ (following the dotted 
arrows). As R is decreased the stable and the unstable 
get closer. Beyond some point {R < Rbc) orbits inside the 
chaotic set do not escape, which means that the chaotic set 
has become an attractor. The right panel shows the attractors 
occurring for R into the coexistence interval. 

giving as a result an (approximately) one dimensional 
curve for the T'^ attractor. We must take an additional 
return map to reduce the T'^ to a fixed point, and the 
chaotic set of dimension (about) four to a line. Consid- 
ering the set of maxima of Zo{t), {Mn}, we took the sub- 
set of maxima whose neighboring maxima were smaller: 
{Nj} = {Mn > Mn±i}, what effectively amounts to con- 
sider the maxima of the low frequency oscillations. The 
results for two values of the parameter R at both sides 
of the crisis (i? = Rbc) are shown in Fig. [9] The chaotic 
attractor exhibits a rough A-shaped structure as occurs 
with the Lorenz map. Probably, the existence of the 
residual fast component makes the attractor to deviate 



significantly from one-dimensionality. In the light of the 
paper by Yorke and Yorke [32i], who studied the transi- 
tion to sustained chaotic behavior in the Lorenz model 
with the Lorenz map, it is found that our results are con- 
sistent with a boundary crisis mediated by an unstable 
three-torus. 



VI. ROUTE TO CHAOS: THEORETICAL 
ANALYSIS 

Condensing all the information obtained from numer- 
ical experiments in the previous sections, we suggest the 
route to high-dimensional chaos represented in Fig. [TUl 
(where it is to be understood that the customary cross 
section through the fast rotating wave is applied). The 
high dimension of our attractor makes somewhat convo- 
luted a geometric visualization. As we explained above, 
we postulate a chaotic attractor whose structure may be 
simplified in terms of a 4-dimensional branched manifold, 
and therefore the Poincare section reduces the attractor 
to a 3-dimensional branched manifold. Figure [TO] repre- 
sents a projection onto K'^, hence some (apparently) for- 
bidden intersections between trajectories appear because 
of the branching. As occurs with the Lorenz attractor 
when it is projected on (say x — z), the intersections 
between trajectories coming from different lobes, and also 
of them with the z axis (that belongs to the stable man- 
ifold of the origin) are unavoidable. Recall that it is the 
moment of the jump when the extra dimension is needed, 
and this is provided by the definition of a branched man- 
ifold. 



Summing up our previous numerical findings, we have 
for descending values of R: the centered PRW (a) be- 
comes unstable through a pitchfork bifurcation (a— >b) 
and two symmetry-related PRWs appear (b). At a su- 
percritical Hopf bifurcation (b— >c) the appears. When 
R is slightly decreased the 2D-torus becomes focus-type 
(d) ; the leading Lyapunov exponent becomes degenerate 
as may be seen in Fig. [2j Hence the unstable manifold of 
the asymmetric PRW forms a 'whirlpool' [i^ when ap- 
proaching the . At Rhet a double heteroclinic connec- 
tion between the asymmetric PRWs and the symmetric 
one occurs (e) (only the upper half of the connection is 
shown for clarity). At this point the chaotic set is created, 
that includes a dense set of unstable 3D-tori. In (f) the 
two simplest 3D-tori are represented with dotted lines, 
because of the heteroclinic birth one of the frequencies of 
these tori is very small. Note that the plot shows that the 
unstable manifold of one of the asymmetric PRWs inter- 
sects the stable manifold of the symmetric PRW which 
in principle violates the theorem of existence and unique- 



ness. This occurs, as we said above, because we are pro- 
jecting the Poincare section onto R'^ (one could imagine 
this as our particular Flatland 61]. Twin secondary Hopf 
bifurcations (f^g) render unstable the 2D-tori and give 
rise to two stable (g). When R is further decreased 
the asymmetric PRWs are not connected by their unsta- 
ble manifolds to the T'^ (h), and the chaotic set becomes 
attracting. 

A two-dimensional cut of the schematic shown in 
Fig. [TO] is depicted in Fig. [11] Only some of the sub- 
plots of both figures are one-to-one related. Specifically, 
Fig. [TlTh) corresponds to R = Rtc, whereas Fig. [TlTi) 
corresponds to i? = R^n ■ 

It is to be stressed that the existence of T'^ attrac- 
tors is not a fundamental part of the transition to high- 
dimensional chaos. Just focus type are needed such 
that the unstable manifolds of the asymmetric PRWs 
form whirlpools [46]. In this way, regarding the chaotic 
attractor, no fundamental change occurred if the unsta- 
ble T'^ shrank to collide with the stable in a subcritical 
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FIG. 10: Three-dimensional representation ol the proposed heterochnic route to create the high-dimensional chaotic attractor. 
Black and gray points correspond to stable and unstable fixed points (cycles in the global phase space), respectively. This 
representation is obtained through a Poincare cross section of the real attractors of the system, that eliminates the (fast) 
oscillation along the ring. 




FIG. 11; Two-dimensional representation of the proposed het- 
erochnic route to create the high- dimensional chaotic attrac- 
tor. It represents a vertical cut of Fig. 1101 (resp. ®) sym- 
bols represent stable (resp. unstable) periodic orbits (2D-tori 
in the full phase space). 

Neimark-S acker bifurcation. This picture would be more 
similar to the transition in the Lorenz system where C± 
become unstable through a subcritical Hopf bifurcation. 

VII. FURTHER REMARKS 

In this section we want to address some subtle aspects 
concerning the transition to chaos shown in this paper. 

The first point that must be noted is that Fig. [10] is 
representing a map (instead of a continuous dynamical 
system) because a Poincare section of the fast dynam- 
ics involving the spatial mode is assumed. However the 



fast spatial frequency is somehow 'orthogonal' along the 
transition and therefore phase space may be understood 
as a direct product of this frequency with the transition 
shown in Fig. [TO] that can be assumed not to differ too 
much from a continuous system. The fast spatial rotat- 
ing wave is approximately conserved along the transition 
(cf . Sec. IIII|) , and the weak interaction of this frequency 
stems from the rotational symmetry of the system and 
manifests through an additional vanishing Lyapunov ex- 
ponent in the high-dimensional chaotic region. 

Let us consider further which is the meaning of the 
Lyapunov spectrum of the chaotic region. According to 
the route described by Figs. [10] and [Tl] an infinite num- 
ber of unstable are created at the heterochnic explo- 
sion (recall that only the two simplest ones are shown). 
Hence, we could expect to have a chaotic attractor with 
three vanishing Lyapunov exponents. Why only two are 
observed? (cf. region CRW in Fig. [2). 

One could be tempted to think that that due to the 
non-hyperbolicity of quasiperiodic behavior, a finite ra- 
tio of the 3D-tori set corresponds to locked tori (prob- 
ably with large denominators), contributing to shift the 
vanishing Lyapunov exponents from zero. However we 
believe that the absence of a third vanishing Lyapunov 
exponent is due to the following. It must be noted that 
the heterochnic connection between the symmetric and 
the asymmetric PRWs is structurally unstable. Nonethe- 
less, we have numerically observed (see Fig. [5]) that the 
connection (approximately) persists. Also, the axial sym- 
metry of Fig.[Tn]has not a theoretical justification. In our 
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view the situation is similar to the study of codimension- 
m (m > 1) bifurcations, where an increase of the highest 
order n < rimax, of the terms constituting the normal 
form (or n-jet) or adding a non symmetric term, makes 
substructures in the bifurcation sets emanating from the 
codimension-TO point to appear. For example, the effect 
of including non-axisymmetric terms to the normal form 
of the saddle-node Hopf (codimension-two) bifurcation 
has been studied by Kirk [47| . 

Hence we expect that, a perturbation on the mech- 
anism shown in Fig. [TU] will break the symmetry that 
allows such a "simple" picture. Inspired on previous 
works [47I l48t dealing with the effect of non-symmetric 
terms on codimension-two bifurcations, we postulate that 
homoclinic connections would replace heteroclinic con- 
nections. A double ('figure-eight') homoclinic of the 
symmetric PRW as well as homoclinic connections of 
the asymmetric PRWs would occur. Symmetric and 
asymmetric PRWs are both saddle-foci so we computed 
their saddle indices finding that a (Shilnikov) chaotic set 
should arise [62] . 

The existence of homoclinic chaos (with the only novel 
feature of having a superimposed fast spatial wave) 
could be considered as uninteresting and well known (see 
e.g. [4^). But it is to be emphasized that as the ex- 
act mechanism is very related to what shown in Fig. 1101 
the first negative Lyapunov exponent is very close to zero 
which makes the information dimension to be larger than 
four (or larger than three if the spatial oscillation is con- 
sidered extra and/or trivial). 

The next question is how the two simplest unstable 
3D-tori appear if the route is not exactly as appears 
in Fig. [ini A possibility is that they appear though a 
saddle-node bifurcation between two unstable 2-tori (see 
footnote 6 for its plausible origin), but we have no way 
to find out this. 

It is well known, [il, [13, HII that global bifurcations 
and complex behavior have in many cases a local ori- 
gin, namely a high codimension point, at which the loci 
of several bifurcations meet in parameter space. In par- 
ticular, a Gavrilov-Guckenheimer (or saddle node-Hopf 
point) yields, among other, quasiperiodic dynamics. The 
refiection symmetry of the model studied in this work 
would imply the possibility of a pitchfork-Hopf inter- 
action in this case. Although the unfolding for this 
codimension-2 point has not been fully characterized, it 
could be a clear candidate as an organizing center of part 
of the complex behavior discussed in this work. How- 
ever, the location of such hypothetical critical point is 
rendered more difficult by the fact that, due to the pres- 
ence of a spatial frequency in all the relevant parameter 
region, the hypothetical pitchfork-Hopf point would oc- 
cur in a Poincare cross section of the system, making all 



the analysis much more involved. In particular, we have 
tried this avenue of research without success, and, thus, 
we have not been able to find a coalescence of the loci of 
the pitchfork and Hopf bifurcations. 

VIII. CONCLUSIONS 

Direct transitions from quasiperiodicity to chaos has 
been observed experimentally (e.g. in [52l|). In the past 
this has been interpreted as an effect of noise or lack of 
good control of the system. Our results indicate that 
high-dimensional chaos may be reached directly without 
the need of noise, but only thanks to a particular type of 
global bifurcation. 

In this paper we have studied by numerical and the- 
oretical arguments the transition to high-dimensional 
chaos, in a system of three coupled Lorenz oscillators. 
The transition from a periodic rotating wave to a chaotic 
rotating wave has been investigated. The structure of the 
global bifurcations between cycles, underlying the cre- 
ation of the chaotic set, is such that a chaotic attractor 
with dimension d w 4 emerges. The transition is not me- 
diated by low-dimensional chaos. Also it must be noted 
that even if the fast rotating wave that is present all along 
the transition is omitted, we still have the creation of an 
attractor with dimension Di > 3. As occurs with the 
Lorenz system, the existence of refiection (Z2) symmetry 
seems to play a fundamental role. 

The high-dimensionality of the chaotic attractor is not 
associated to hyperchaos. Far from it, there is only one 
positive Lyapunov exponents but high-dimensionality is 
possible due to the existence of two vanishing and one 
slightly negative LEs. Hence, according to the Kaplan- 
Yorke conjecture the information dimension is above 
four. We have also measured the correlation dimension 
obtaining a value very close to four. The degeneracy of 
the null LE make us think that a set of unstable tori is 
embedded in the attractor. 

We have focused on giving a geometric view of the bi- 
furcations occurring in the 9-dimensional phase space of 
the system. Although the precise seq uence of bifurcations 
is probably resistant to analysis [63j, we have been able 
to give a geometric view of the transitions that explains 
the emergence of the chaotic set, through a 'heteroclinic 
explosion', and its conversion in attractor. This step oc- 
curs through a boundary crisis when the chaotic attractor 
collides with its basin boundary formed by two unstable 
3D-tori. In consequence a power law for the mean length 
of the chaotic transients is observed. 

This work was supported by MICINN (Spain) and 
FEDER (EU) under Grants No. BFM2003-07749-C05- 
03 and FIS2006-12253-C06-04 (DP), and FIS2007-60327 
(FISICOS) (MAM). 
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